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INTRODUCTION 


The  need  for  more  advanced  optimal  control  techniques 
has  become  apparent  in  recent  years.  With  the  advent  and 
perfection  of  minicomputers  (microprocessors)  digital 
aircraft  control  has  come  of  age.  Need  for  control 
optimization  can  be  applied  to  commercial  aircraft  in  the 
form  of  fuel  economizing  and  increased  desire  for  instrument 
landings,  as  well  as  its  current  application  to  defense  and 
aerospace  needs.  One  major  advantage  of  this  control 
technique  is  the  reduction  or  elimination  of  the  need  for  a 
pilot,  for  certain  tasks,  thus  reducing  the  possibility  of 
human  errors  while  retaining  the  effects  of  human  judgement. 

In  order  to  apply  optimal  control  techniques,  certain 
assumptions  and  simplifications  are  usually  made.  One  such 
simplification  involves  the  linearization  of  the  plant 
equations  around  one  flight  condition.  The  question  then 
arises  as  to  how  this  will  affect  the  response  of  the  system 
at  other  flight  conditions.  Using  sensitivity  analysis  this 
thesis  examines  these  effects  along  with  the  effects  of 
using  an  observer  to  estimate  certain  unavailable  states. 
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The  method  used  here,  begins  by  deriving  the  equations 
of  motion  of  an  airplane.  These  equations  are  used  to 
simulate  aircraft  flight.  Using  this  model  as  the  plant  to 
be  controlled,  the  nonlinear  equations  are  linearized  around 
a soecified  flight  condition.  With  these  linearized 
matrices  a PID  control  is  then  designed  by  using  the  minimum 
principle.  For  the  case  where  certain  states  are  not 
available,  a low  order  observer  is  then  designed.  This 
controller  can  then  be  applied  to  the  nonlinear  system  and 
tested  for  its  acceptability. 

The  objective  of  this  thesis,  therefore,  is  to  study 
the  practicality  of  this  control  under  simulated  flight 
conditions.  This  is  accomplished  first  by  determining  how 
sensitive  the  closed  loop  system  is  to  various  parameter 
variations.  Then  the  effects  of  using  an  observer  are 
examined.  Finally,  the  controls  are  tested  to  see  how  well 
they  perform  in  an  actual  real  time  application. 

The  procedure  followed  in  this  thesis  is  typical  of  the 
method  used  for  design  of  a linear  control  for  application 
to  a nonlinear  plant.  The  derivations  are  general  enough 
that  this  same  procedure  could  be  used  to  control  a variety 
of  other  systems. 
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CHAPTER  2 


AIRCRAFT  MODELLING 


If 


r 


In  order  to  proceed  with  precise  and  efficient 
controlling  of  a system,  an  accurate  model  of  the  system 
must  first  be  found.  In  this  thesis  a simplified  model  of 
an  airplane's  vertical  equations  of  motion  is  used. 
Theoretically,  though,  any  system  could  be  used  and 
controlled  using  the  techniques  of  the  following  chapters 
(e.g.  lateral  equations  of  motion). 


2.1  EQUATIONS  OF  MOTION 


The  equations  necessary  for  this  model  can  be  derived 
using  simple  dynamics  while  making  certain  reasonable 
assumptions  about  the  aircraft  (ignoring  aeroelasticity , 
etc.).  In  general  an  airplane  can  be  assumed  to  have  the 
configuration  as  shown  in  Fig.1  where  the  symbols  refer  to 
the  quantities  as  given  in  Table  I.CRef.  4,6,7). 


The  X-axis  is  usually 
angle  of  attack  (a)  is 
approximations  can  be  made, 
sin (a  ) = a 
cos (a  ) = 1 


(by  definition)  such  that  the 
small.  Therefore  small  angle 
This  leads  to  the  following. 


u = Vcos (a  ) =V 


Figure  1.  Airplane  Configuration 

Table  1 

Symbols 

a:  angle  of  attack 
9 : pitch  angle 
y:  flight  path  angle 
M:  mass  of  the  aircraft 
V:  veloci-ty 
H:  altitude 

W:  weight  of  the  aircraft 
I yy : moment  of  inertia 
XCg : center  of  gravity 
S:  wing  surface  area 
p : air  density 
c:  chord  length 
X:  body  axis 

Z:  vertical  normal  to  body  axis 


■1.  J » Ik  ■*!... UM. WJ 


u=Vcos(a)-Vsin(a)a=V-Vaa 

w=Vsin(a)=va 

w=Vsin(a)+Vcos(a)a=Va+Va=Va 

sin(Y)=sin(Y-a)=sin(V)cos(a)-cos(y)sin(a) 

= sin ( Y ) -acos (Y  ) 

Summing  the  forces  in  the  X-direction : 

0=-MV+T-D-Wsin(Y)+Lsin(a)-(D-T)cos(a) 

Now  using  these  approximations: 

0 = -MV  + T - D - W(  sin(Y)  - acos(Y)  ) 

+ (L  - Wcos(Y)  - MV(0-a)  ) 

So,  using  2.2)  from  below: 

2.1)  0 s -MV  + T - D - Wsin(Y ) 

Summing  the  forces  in  the  Z direction: 

0 = -M(w  - uq ) + Wcos(Y)  - Lcos(a)  - (D  - T)sin(a) 


Drag  and  thrust  terms  are  generally  smaller  than  the 
lift  and  weight  terms,  and  about  equal  to  each  other.  Using 
this  plus  the  fact  that  they  are  multiplied  by  a small 
angle,  a,  the  last  term  above  is  considered  negligible.  So 
this  simplifies  to: 

2.2)  0 = L - Wcos(Y)  - MV(0-a) 

Summing  the  moment  in  the  X-Z  plane: 

2.3)  19  = My 


Summarizing  these  equations: 

2.4)  V = -jJ-CT  - D - Wsin(Y)  ) 


Q - yfV  (L  - Wcos  (Y  ) ) 


a = 


6 


2.5)  H = Vsin(y) 

Wy/I>y 

L = pj  V25C1 
D = ^pV2SCd 
My  = jOV2S<=Cm 

Where  Cj,Cm,  and  depend  on  wing  plan  form  used  and 

placement  of  the  wing  (and  sometimes  placement  of  the 
engines).  All  the  coefficients  in  these  equations  can  be 
found  for  any  size  airplane  using  the  specified 
configuration  and  by  looking  up  the  wing  specifications  in  a 
book  listing  wing  characteristics.  These  equations  are 
generally  simplified  for  Mach  numbers  less  than  1.0  by: 


C1 

= Clo 

+ 

Cla  + Clf8f 

cd 

= Cdo 

+ 

Cdcl  + Cdf8f 

Cm 

= ^mo 

+ 

^mcl^l  + ^me 

a + q) 

Any  airplane  now  can  theoretically  be  simulated, 
perhaps  with  minor  modifications  due  to  engine  placement, 
tail  configuration  or  Mach  number.  For  simplicity  the 
coefficients  of  the  GAT  II  simulator  as  described  in  Daly's 
thesis  (Ref.  4)  are  used  with  minor  revisions. 

Thrust  is  a more  complicated  subject.  It  is  highly 
dependent  on  Mach  number  , altitude,  and  the  type  of  engine 
used  (turboprop,  turbofan,  propeller,  etc.).  In  general 
there  are  no  easily  found  formulas  for  thrust.  Again  for 


7 


simplicity,  and  due  to  a limited  knowledge  of  this  subject, 
the  thrust  formulation  (propeller)  used  in  Daly's  thesis  was 
adopted,  which  is: 


Map 

= CPO 

+ CphH 

+ cpnM  + CpntNSt 

Bhp 

II 

O 

c r 

O 

+ cbnM 

+ cbpMa?  + cbhH 

T = 

Me3hp ( 

o 

+ 

o 

c_> 

tvV  + cthH  + CtvhVH) 

where  the  abbreviations  are  as  listed  in  Table  2. 


Table  2 

Thrust  abbreviations 
Map:  manifold  pressure 
N : RPM 

Bhp:  brake  horsepower 
T:  thrust 
5t : throttle 
Me:  number  of  engines 


Throughout  this  work  just  one  flight  condition  was 
examined,  although  any  reasonable  conditions  could  have  been 
examined.  The  values  that  are  used  are  shown  in  Table  3. 

Combining  all  these  equations  yields  the  fifth  order 
nonlinear  system  be  low,  where  the  states  x i - x^  and  controls 
u i -u ^ represent  the  quantities  listed  in  Table  4. 


= x. 


“MxpS  Px|s( 


ulo 


1 a x 1 


:lfu2 ) 


- Woos  (x-j  - x 


— m m 


4 


— TUI  ■ 


Coefficients  and 

flight  conditions  use'd 

ClQ  = 0.0765 

Cp0  = 29.92 

C la  = 4.62 

Cph  = 0.0009 

Clf  = 0.365 

Cpn  = 0.00076 

Cd0  = 0.026 

Cpnt  = °-0l65 

Cdcl  = °-062 

Cb0  = -352.3 

Cdf  = 0.021 

Cbn  = 0.1155 

Cmo  = O'1 

Cbp  = 10.8 

Cmcl  = -0-0529  + £eg 

Cbh  = 0.0025 

Cme  = -0.0354 

Cto  = 3.5 

Cmf  = -0.0368 

Ctv  = -0.00642 

Cbhp  = 2.11 

cth  = -4.73(10)~5 
ctvh  = 8.7(10)“® 

N = 2500  rpm 

VQ  = 190.66  ft/sec 

P = 0.004842  slugs/ f 1 3 

H0  - 2000.  ft 

S = 180.  ft2 

0o  = % = <*c  = o . 

xcg  = 0.2  ft 

c = 5 ft 

W = ^000.  lbs 

Iyy  = 2050.  slugs  ft2 


2.8) 

2.9) 

2.10) 


1 

I 


x2 

= lV(- 

W sin ( x^ 

- 

xi)  - 

7Px2s(cdo  + ■'del 

(C 

lo 

+ C la  x 1 

+ 

Clfu2 

)2 

+ Cdfu2  + (c 

to 

+ Ctvx2 

+ 

Cthx5 

+ 

^tvhx2x5  ^ ^bo 

+ 

cbnN 

+ Cbhx5 

+ 

Cbp<C 

po 

+ CpnM  + Cphx5 

+ 

2 p n t 

x3 

= x 4 

x4 

= ~b 

yy 

(■^“px2Sc ) 

r r 
L t 

mo  + 

^mc 1 ^ C lo 

+ Cla  x 1 

+ 

Clfu2 

) + 

Cmf u2  “2^2(x1 

+ 

x4 ) 

= x2sin(x2  - x 1 ) 


Table  4 

States  and  controls 
x 1 = a u-,  = 5e 

x2  = V u 2 - 8 f 

x3  = 9 u3  = St 

x4  = 6 

xc  = H 


L 


2.2  LINEARIZATION 


In  order  to  use  linear  control  theory, it  is  necessary 


to  obtain  a linearized  version  of  the  system, 


.s  is 


accomplished  using  first  order  perturbation  techniques 
x = f ( x , u , t ) 


if  if 

x = ax  x + au  u 


if 
a x > 


Define : 


3 = 11 
au 


2.3  EQUILIBRIUM 


To  proceed  we  must  now  find  an  equilibrium  point. 
Given  a desired  set  point,  subject  to  the  restrictions 
below,  one  can  find  what  the  control  values  should  be  at 


equilibrium. 


0 = f(X,U) 


where  X and  U are  equilibrium  values. 

Given  the  desired  equilibrium  state  X, 
since  x^  = 0 and  = 0 we  get, 

2.11)  X4  = 0 

2.12)  X3  = X1 


From  x 


0 we  see  that: 


2 • 1 3 ) u2  <-cl0  - claXl  + c|cos(X3-X1)  + mx4x2) 


= rr-.C-C1^  - ClaX1  +g^) 


=q^-^io 


2.14)  Ut  = cr- (Cmclc1  + cni0  + Cmfu2) 


2.15)  U 


3 CbpCpnt^  ^bo  + cbDCpo  + ^^bn  + ^bpCpn^ 
+ (Cbh  + cbpcph)x5  “ T§npCd 


where , 


« = jpxl 

C1  = Clo  + ClaX 1 + clfu2 

Cd  = cdo  + c d c 1 c 1 + cdfU2 

Tbhp=  Ne(Cto  + "tvX2  + CthX5  + CtvhX2X5) 


Pluggirg  in  these  eauilibrium  values  and  linearizing 
yields  the  coefficients  of  A and  B given  in  Tables  5 and  6. 
\'ow  our  linearized  plart  with  feedback  can  be  modeled  as  in 
the  block  diagram  given  in  Fig. 2. 


Table  5 


n 


i = — ■ (C — C 

zyy 


me  I1"  la  “ 2X^a  1 1 ^ 


*42 


= !h^[  (cmo  * Cmc,cl0)  * Cm0lclax,  * cmu, 


yy 


+ (CmclClf  + Crnf)U2 


I 


I 


Table  5 

Linearized  '3'  coefficients 


V> 

Qi  r 

m m 

b12 

MX/1? 

• * 

b22 

= 

-3^2CdclCli 

~c  i + 

C d f ^ 

b23 

s 

ThhpCbpCpnt? 

i 

b4 1 

- 

QScr 
t ume 

yy 

b42 

= 

7 — umcr-lf 
yy 

+ Cmf 

- — b i p ) 
2X2  12 

b11 

= 

b 1 3 = b2 1 = 

b3 1 = 

b32  = b33 

= b 

43 

= b51  = b52 

= b53 

= 0 

U 
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CHAPTER  3 


CONTROLLER  DESIGN 


3. 1 FEEDBACK  CONTROL 

Since  a pilot  is  generally  interested  in  obtaining  a 
desired  velocity,  altitude,  or  pitch  angle  it  would  be 
favorable  to  have  controls  which  directly  relate  to  these 
quantities.  As  it  is,  he  must  estimate  how  much  more 
throttle  he  must  use  to  get  the  airplane  to  a certain 
velocity,  and  similar  situations  exist  for  altitude  control. 
This  is  not  a critical  problem,  usually,  since  there  is 
little  need  for  perfect  accuracy  when  one  is  20,000  ft.  in 
the  air.  But  in  order  to  conserve  fuel  or  more  importantly, 
when  landing  in  a heavy  storm  this  can  become  critical.  One 
other  application  is  for  pilotless  control  of  an  aircraft, 
as  is  the  case  for  missiles  and  drones. 

Assuming  that  the  state  variable  vector, X,  is  available 
(or  at  least  estimated)  a logical  choice  for  a control  is 
one  which  the  control  variables  are  a linear  combination  of 
the  state  variables.  In  other  words  the  control  can  be 
written : 


u = Fx 
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Mow  in  order  to  analyze  the  stability  of  the  system 
with  this  control,  the  following  linearized  system  is 
examined : 

x = Ax  + Bu  =Ax  +BFx=(A  + BF ) x 

3.1.1  Pole  Placement  Technique 

Using  pole  placement  techniques,  the  eigenvalues  of 
this  system  can  be  placed  arbitrarily  if  the  system  is 
controllable  (Ref. 1). In  other  words  if  the  matrix 
( B , AB , A^B , . . . , An“^ B ) has  rank  n,  a desired  response  can  be 
obtained.  Both  the  controllability  test  and  the  placement 
of  eigenvalues  can  be  accomplished  using  Linsys  (Ref. 2). 
The  poles  are  placed  in  the  left  half  plane  so  that  as  t 
gets  large  x=0.  This  character istic  will  be  utilized 
later . 

3.1.2  Optimal  Control  - the  Minimum  Principle 

An  alternate  method  of  finding  the  control  feedback 
matrix,  is  to  assume  a performance  index  J,  that  we  want  to 
minimize.  If  we  assume  our  performance  index  to  be  of  the 
f orm: 

J =j  (xtQx  + uTRu  ) dt 

where  Q and  R are  cost  matrices  which  can  be  adjusted  to 
procure  a certain  response.  If  A,  B,  and  Q conform  to  the 
conditions  stated  in  Appendix  5 then  there  exists  a K,  such 
that  the  system  with  the  following  control  is  asymototica  lly 
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stable  (Ref.  1 ),  where  K is  found  from  the  algebraic  Riccati 
equation: 

-KA  - ATK  + KBR_1BTK  - Q = 0 
and  u = Fx  = -R"1BTKx 

The  above  Riccati  equation  and  the  resulting  F can  all 
be  found  using  Linsys  (Ref.  2 ). 

3.1.3  Augmentation 

Since  both  of  these  above  methods  result  in  systems 
which  are  asymptotically  stable,  we  can  use  this  fact  to 
'force'  a state  or  states  to  a specified  set  point  value. 
We  note  that: 

x (°° ) = 0 
y = Cx 

If  we  now  want  y(»)  = yref  = Cx,  we  can  augment  our 

original  system  and  define  x'  such  that: 


x 

= 

’A 

d 

x ’ 

B 

u 

0 

y-yref 

C 

°J 

0 

. _yref 

or 

x'  = A'x'  + B'u  + v 
and  with  proper  feedback: 

x'(«)  = 0 and  y(°°)  = yref 

In  our  case  we  want  to  control  velocity,  V,  and 
pitch  angle , 0 , so : 

*6  = x2  “ vref 
*7  = x3  " *ref 


3.2  OBSERVER  DESIGN 


Throughout  the  rest  of  this  thesis  we  will  assume  that 
certain  states  are  not  actually  available  for  measurement 
(or  perhaps  not  very  accurate).  Since  for  accurate  and 
efficient  controlling  (or  perhaps  for  pilots  reference) 
these  states  are  needed,  some  technique  is  necessary  to 
provide  these.  For  this  an  'observer'  or  'state  estimator' 
is  used  which,  as  the  names  imply,  estimate  the  state  or 
states  that  can't  be  directly  measured  . For  the  purposes 
of  this  thesis  (and  to  demonstrate  this  technique)  it  will 
be  assumed  that  the  states  that  are  not  available  are  angle 
of  attack  and  pitch  rate.  In  reality  these  may  be 
measurable,  but  for  reasons  such  as  cost  or  precision  it 
will  be  assumed  that  they  are  not  accurately  available.  On 
the  other  hand  the  3 states  - pitch  angle,  velocity,  and 
altitude  have  to  be  measured  for  the  pilots  referance, 
anyway,  so  these  certainly  would  be  available,  and  would  not 
require  any  further  instrumentation. 

There  may  be  a realistic  limit  on  the  order  of  the 
observer,  due  to  computer  size,  speed  or  cost,  so  a low 
order  observer  will  be  used.  Using  the  results  obtained  in 
Appendix  A: 

l1 1 
L21 
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w2 

- G 1 W2  + G2y 

+ G^u 

1 • 

where : 

<*1 

= a22  - ka12 

1 

■« 

G2 

= A21  - KAn 

+ G-,K 

g3  = 

b2  - 

KB1 

w2  = 

A 

z2  - 

Kz  1 

Because 

of 

the  special  form  of 

the 

transf  ormation 

matrix,  T, 

the 

actual  implementation 

is 

made  somewhat 

simpler.  This  is  partly  due  to  the  fact  that  3 of  the 
states  are  outputs,  ar.d  part  due  to  the  choice  of  M.  The 


matrix,  T,  is 

chosen 

as : 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

T = 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

As  is  easily  seen  x is  obtained  simply  by  interchanging 
the  order  of  z^.  So  the  estimated  states  z 2 are  actually  x^ 


and  x2,  respectively. 

It 

might 

also  be  noted 

here  that  it 

is  assumed  that 

x6 

and 

Xj  are 

available  since  they 

are 

needed  previously 

f or 

the 

augmentation.  Taking  this 

into 

account  results 

in 

the 

block 

diagram  of 

F i g . 3 . 

The 

nonlinear  version 

is 

shown 

in  Fig. 4, 

where 

the 

observer-controller  is  the  same  as  in  the  linear  case, 
except  that  the  controls  and  states  are  normalized  around 

f * 

ii 


the  equilibrium  point . A slightly  different  approach  can 
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P 

I • 
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2.  Relatively  smooth  response,  with  few,  if  any,  jumps 
°r  Jerks' 

3.  As  little  overshoot  as  possible. 

4.  Keeping  controls  within  bounds  and  limitations. 
Following  this  set  of  conditions  the  cost  matrices  were 


arrived 

at . 

Upon  substitution 

into  the 

Riccati  package  of 

Linsys 

the 

feedback  matrix, 

F , was 

calculated.  The 

eigenvalues  for  this  linear  system,  along  with  the  rest  of 
the  system  matrices,  are  shown  in  Appendix  D.  It  should  be 
noted  here  that  in  order  to  use  the  Riccati  package  the 
matrices  were  normalized  so  that  all  elements  were  about  of 
the  order  of  one  (X2»  Xg,  xg  are  scaled  down  by  100). 

The  observer  feedback  matrix  was  similarly  chosen  using 
the  Riccati  package.  The  values  obtained  for  K and  the 
corresponding  Q,  R,  and  eigenvalues  are  given  in  Appendix  D. 


3.3.2  Effects  of  Parameter  Variations 

In  order  to  test  how  effective  the  controls  are  in 
simulated  flight  conditions,  several  types  of  parameter 
variations  were  used.  All  of  the  variations  are  realistic 
disturbances  that  an  aircraft  can  encounter.  The  time 
responses  due  to  these  variations  appear  in  Figs.  5. 1-5.7  on 
pages  24-34. 
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3.3.2a  Perturbations 

The  first  type  of  disturbance  studied  was  simply  a 
change  in  initial  condition,  which  could  represent  a gust  of 
wind,  for  instance.  As  one  can  see  from  the  plots  of  Figs. 
5.1a  and  5.2a,  the  controller  (with  observer)  performed 
properly  as  long  as  these  perturbations  were  within 
reason(i.e.  within  30  ft/sec  or  .2  rad). 

3.3.2b  Disturbance 

For  this  a constant  disturbance  vector  is  introduced 
into  the  plant's  differential  equations,  this  could 
correspond  to  a headwind  or  perhaps  inaccurate  readings. 
Figs.  5.3a  and  5.3c  show  the  responses  with  this  disturbance 
introduced.  For  a small  disturbance  the  response  was  hardly 


af f ected 

. As 

the 

disturbance 

was  increased  the 

response 

showed 

signs 

of 

instability 

and  the  response 

became 

oscillatory. 

3.3.2c  Set  Point 

The  set  point  is  the  desired  flight  condition  as 
specified  by  the  pilot  (e.g.  velocity  or  pitch  angle).  An 
example  of  this  is  shown  in  Figs.  5.4a  and  5.4c.  For  set 
points  far  from  the  original  conditions  the  response  became 
unfavorable  and  often  the  controls  tended  to  exceed  their 
realistic  limits. 
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3. 3. 2d  Plant  Variations 

It  was  also  possible  to  change  the  different  aircraft 
parameters,  such  as  center  of  gravity,  to  represent  a shift 
of  cargo  as  in  Figs.  5.6a  and  5.6b,  for  instance.  Bringing 
the  center  of  gravity  farther  forward  had  little  effect 
until  around  .5  ft.,  where  depending  on  whether  an  observer 
was  used  or  not  the  system  went  unstable,  with  the  observer 
system  going  unstable  earlier.  Bringing  the  center  of 
gravity  farther  back  in  the  plane  tended  to  yield  better 
responses  than  at  the  nominal  condition,  and  evidently 
created  a more  stable  situation.  This  agrees  with  what  is 
expected  from  aircraft  dynamics,  the  farther  back  the  center 
of  gravity  the  more  stable  it  should  be  by  design. 

3.3.3  Effects  of  Using  an  Observer 

The  use  of  the  observer  did  not  seem  to  drastically 
affect  the  response  of  the  system.  As  can  be  seen  from  the 
plots  of  Fig. 5.1  there  is  no  significant  difference  after 
the  first  2 seconds.  The  observer  system  did  seem  to  go 
unstable  sooner,  as  variations  approached  their  reasonable 
limits.  But  the  observer  also  tended  to  'smooth'  the 
response  somewhat,  which  is  a desirable  side  effect  of  its 
use . 
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3.3.4  Real  Time  Results 

To  simulate  actual  aircraft  implementation,  the 
controller  designed  in  this  chapter  was  tested  on  the 
AD5-PDP11  hybrid  computer.  The  continuous  model  was  used 
here  rather  than  convert  to  a discrete  model.  This  is  seen 
to  be  a safe  simplification  for  high  enough  sampling  rates. 
Examples  of  the  responses  are  shown  in  Figs.  5. 1,5. 3,  and 

5.4  ( b and  d ) . 

3.3.4a  Sampling  Rate  Effects 

The  system  seemed  to  be  very  sensitive  to  the  sampling 
rate.  It  was  found  that  the  digital  controlling  program 
took  around  .01  seconds  to  run.  But  the  system  with 
controller  tended  to  go  unstable  with  .08  second  sample 
time.  This  is  a rather  limited  range  of  sampling  rates  and 
would  tend  to  be  a drawback  for  real  time  operation. 

3.3.4b  Observer  Effects 

Here,  as  was  the  case  for  the  continuous  time 
simulation,  the  observer  use  had  little  effect  on  the 
response.  It  tended  to  also  smooth  the  response  somewhat 
which  helped  to  keep  the  controls  in  bounds.  It  was  also 
noted  that  the  observed  states  had  little  effect  in  the 
first  place.  This  was  discovered  by  not  feeding  these 
states  back  at  all. 


i sure  5.2a  Stats  Response  with  Observer 
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3.3.4c  Real  Time  Effects 

The  responses  were  somewhat  more  jumpy  probably  because 
of  the  analog  to  digital  conversion  which  would  reduce  the 
accuracy,  especially  in  the  case  of  altitude.  In  general, 
though,  the  responses  in  real  time  were  nearly  identical  to 
the  continuous  time  counterpart,  and  if  anything  tended  to 
have  less  overshoot. 

3.3.5  Operating  Envelope 

When  conditions  were  varied  too  far  from  the  nominal 
flight  conditions  the  response  character istics  deteriorated. 
With  this  in  mind,  different  pertinent  parameters  were 
varied  to  determine  a region  for  which  this  controller  would 
still  operated  satisfactorily.  This  operating  envelope  is 
found  to  be  approximated  by  the  following  boundaries  with 
the  corresponding  characteristics  making  the  response 
undesirable . 

1.  Velocity  170<V<240.  At  170  ft/sec  the  pitch  angle  is 
oscillatory.  At  240  ft/sec  the  throttle  control  goes  out  of 
bounds . 

2.  Pitch  Angle  -0.25<a<0.1.  The  throttle  goes  out  of 
bounds  at  both  of  these  limits. 

3.  Disturbance,  velocity:  10  ft/sec^  pitch  angle:  0.2 
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rad/sec.  The  throttle  goes  out  of  bounds  for  both  of  these. 

4.  Center  of  Gravity  -3 . °<Xcg< . 25 . At  Xcg  = 0.25  the 
system  becomes  oscillatory.  At  = -3.0  the  response  is 
sluggish  and  the  controls  go  out  of  bounds. 

5.  Air  Density  0 . 003< P<0 . 0 1 . At  the  lower  bound  the  system 
oscillates  while  at  the  upper  limit  the  throttle  goes  out  of 
bounds . 
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CHAPTER  4 


SENSITIVITY 


Sensitivity  is  the  analysis  of  how  a parameter 
variation  affects  the  dynamics  of  a system.  For  instance, 
by  shifting  an  aircraft's  center  of  gravity  (e.g.  fuel  or 
cargo  moving)  too  far  forward,  an  inherently  stable  aircraft 
can  become  unstable,  making  it  harder  to  control.  So  one 
purpose  of  feedback  control  is  to  hopefully  reduce  or 
eliminate  any  tendency  toward  instability.  In  this  chapter, 
the  sensitivity  of  the  controlled  system  and  the  observer 
controlled  system  are  compared. 

4.1  SENSITIVITY  MODEL 

By  definition,  trajectory  sensitivity  is  the  change  in 
trajectory  due  to  a change  in  the  parameter  p. 

= & 

Or  normalizing  this  with  respect  to  p: 

lim  X§  = p 3x 
Sp 

If  we  are  given  a system  of  the  form: 
x = Ax  + Eu 

taking  the  partial  derivative  of  this  eauation  with  respect 
to  p yields: 


Here  we  assume  that 


u does  not  deoer.d  on  o. 


u 


utilizes  feedback  control,  u = Fx,  then  we  assume  A' 
BF ) , where  F does  not  depend  on  p.  Letting 
4.1  ) 


= (: 


JX  _ A > dx  . b\  A3,, 
^p-A5p  + dpX  + 3pU 


i ' dx 

dD 


( AA  , A Btr ) ,, 

C§p  + ¥pL  ;x 


Now  we  have  an  n^'a  order  sensitivity  inde  x , p , which 
requires  x and  u as  its  inputs  or  for  simple  feedback 
control  just  requires  x.  For  this  system,  the  block  diagram 
is  as  giver  in  Fig. 6. 


Figure  6.  .System  Sensitivity  Block  Diagram 


If  an  observer  is  used  this  complicates  the  situation 
considerably,  because  instead  of  direct  feedback  control  we 
use  the  observed  value's  feedback: 
u = F x = FT  z 

We  assume  that  F and  7 do  not  depend  directly  on  o,  so: 


— 
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The  job  now  is  to  find  how  z depends  on  p, 
observer  chapter  we  know: 


From  the 


A 

Z = 


z2 


$2  = W2  + Kz 


As  before  |§1  is  available  simply  in  the  form  of  |~. 


dV 


* ap> 

U2  = 1^2  + K 


ap 


3P 


ap 


a v 
k a p 


So 


aw 


2 is  now  needed. 


W2  = G -j  W2  + G2y  + G^u 

So  taking  the  partial  derivative  with  respect  to  p: 

4.2)  |g2 

which  is  an  (n-p)Ln  order  system  (requiring  integration). 
Converting  this  into  a block  diagram  yields  the  system  of 
Fig. 7. 


Gl|S2  * g2SF 

, th 


r ^ u , a G 1 
G33P  +&p1w2 
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4.2  SENSITIVITY  RESULTS  * 

Before  proceeding  with  examining  the  sensitivity  plots 
it  is  important  to  check  the  curves,  since  they  are  useless 
unless  they  are  giving  the  correct  information.  Also  this 
gives  a better  understanding  of  their  meaning.  Fig.  8.6 
represents  a negative  shift  of  25%  of  the  center  of  gravity 
from  that  in  Fig.  8.1a.  Concentrating,  for  instance  on  the 
overshoot  peak  of  pitch  angle,  at  about  1.2  sec.  there  is  a 

* Figs.  8. 1-8.8  appear  on  pages  ^4-52  where  0 refers  to 

, while  R refers  to  f-£. 
eg 
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difference  of  about  -.004.  Mow  turning  our  attention  to  the 
corresponding  sensitivity  curve  (Fig.  3.1b)  of  pitch  angle 
at  1.2  sec.  the  value  is  around  0.016.  When  multiplied  by 
the  negative  25*  shift  yields  -0.004  which  checks  with  what 
was  obtained  above.  The  rest  of  the  curves  can  similarly  be 
checked  in  this  manner  for  accuracy.  The  remainder  of  this 
chapter  deals  with  interpreting  the  results  of  the 
sensitivity  analysis. 


4.2,1 


Perturbation  Effects 


From  Figs.  8.1a  and  8.1b  or  Figs.  8.2a  and  8.2b  the 
trajectory  for  the  pitch  angle  sensitivity  is  in  phase  with 
the  trajectory  for  pitch  angle  itself.  This  can  be  seen  by 
noting  for  instance  that  during  approximately  the  first  1.5 
seconds  the  sensitivity  is  positive  while  the  pitch  angle  is 
also  primarily  positive  ignoring  the  initial  transition 
period.  In  other  'words  if  the  center  of  gravity  is  shifted 
forward  (in  the  positive  direction)  the  effects  of 
perturbation  will  be  accentuated.  After  the  first  3 seconds 
some  oscillation  is  also  introduced.  Taking  this  into 
account  it  can  be  deduced  that  a positive  shift  of  center  of 
gravity  will  create  a less  stable  situation.  Cn  the  other 
hand  velocity  seems  to  yield  a more  favorable  response  for 
the  same  condition.  This  is  concluded  by  noting  that  during 
rise  time  (the  first  2.5  seconds)  the  trajectories  for 
velocity  and  its  corresponding  sensitivity  are  in  phase, 
while  during  the  overshoot  period  the  phases  are  opposite. 


i 


42 


it 


Thus  a quicker  rise  time  and  less  overshoot. 

4.2.2  Disturbance  Effects 

From  Figs.  8.3a  and  8.3b  the  primary  effect  of 
shifting  the  center  of  gravity  forward  is  to  increase  the 
overshoot  for  velocity  but  decrease  it  for  pitch  angle  as  is 
indicated  by  the  plots  using  similar  reasoning  as  above. 

4.2.3  Set  Point  Effects 

Examining  Figs.  8.4a  and  8.4b  shows  that  a shift  of 
center  has  very  little  effect  on  the  response  with  respect 
to  set  point  specification.  This  is  determined  by  the 
relatively  small  magnitudes  of  the  sensitivity  trajectories. 
What  little  effect  there  is  gives  the  pitch  angle  a quicker 
rise  time  (sensitivity  positive  in  the  first  2 seconds)  and 
reduces  the  overshoot  somewhat  after  that.  There  is  no 
specific  effect  noticed  on  velocity. 

4.2.4  Equilibrium  Point  Effects 

The  sensitivity  does  not  seem  to  be  significantly 
affected  by  operation  at  an  equilibrium  point  other  than  the 
nominal,  as  is  indicated  by  comparing  the  sensitivities  of 
Figs.  8.1a  and  8.6. 
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4.2.5  Observer  Effects 

The  observer  appears  to  have  favorable  effects  on  the 
response  of  both  velocity  and  pitch  angle.  As  can  be  seen 


f rom  Figs.  8.1a 

and 

8.1c, 

both 

velocity 

and  pitch 

angle 

sensitivity  are 

in 

phase 

with 

their  cor 

responding 

state 

variable  during 

rise 

time  , 

whi  1 

e in  the 

overshoot 

period 

they  are  of  opposite  phase.  Thus  the  effect  of  the  observer 
is  to  speed  up  rise  time  while  reducing  overshoot  somewhat. 
The  observer  also  seems  to  be  somewhat  more  sensitive  than 
without  an  observer,  since  the  curves  are  of  considerably 

| ! 

larger  magnitude  than  those  without. 

4.2.6  Air  Density  Sensitivity 

A 

. * Examining  Figs.  8.1a  and  8.8a  or  8.2a  ar.d  8.3b 

ind- cates  that  for  a higher  (positive  shift)  air  density  the 
response  to  perturbation  gets  reduced.  This  is  because  when 
the  state  is  positive  the  corresponding  sensitivity  is 
negative,  tending  to  reduce  the  effect  of  the  perturbation. 
The  effects  of  disturbance  inputs  are  different.  From  Figs. 
8.3a  and  3.8c  the  overshoot  periods  for  both  velocity  and 
pitch  angle  seem  to  get  lengthened  (sensitivities  are  in 
phase  and  extend  longer)  for  a positive  shift  of  air 
density.  For  pitch  angle  the  initial  the  initial  overshoot 
magnitude  gets  decreased  while  for  velocity  it  gets 
increased.  Comparing  the  sensitivities  with  and  without  the 
• observer  (Figs. 8. la  and  8. Id)  the  observer  amplifies  the 

effects  of  changing  air  density. 


System  Sensitivity  with  Observe 


In  this  thesis  just  a few  aspects  of  using  optimal 
control  theory  as  applied  to  aircraft  flight  have  been 
examined.  From  this  many  desirable  characteristics  were 
revealed.  But,  on  the  other  hand  there  were  several 
shortcomings,  which  would  have  to  be  resolved  before  this 
method  can  be  fully  accepted  and  utilized. 

I 

A review  of  this  control  method  shows  some  valuable 
features  that  should  not  be  discounted.  For  example, 
controlling  of  the  aircraft  involves  only  simple  command 
responses  as  opposed  to  a need  for  good  knowledge  of  flight 
controlling,  so  pilot  error  can  be  greatly  reduced  or 
eliminated.  Along  these  same  lines  near  perfect  accuracy 
can  be  attained.  Furthermore  a desired  response  can  be 
chosen  such  as  for  maneuvarabi lity  or  smoothness  of  flight, 
with  automatic  disturbance  rejection. 

Most  of  the  drawbacks  encountered  had  to  do  with  flying 
outside  of  the  nominal  flight  conditions.  For  instance,  the 
plane  would  go  unstable  for  a relatively  small  forward  shift 
of  center  of  gravity.  This  problem  could  be  reduced  by 
originally  designing  the  controls  around  a more  forward 
center  of  gravity,  since  the  system  stayed  stable  for  larce 


I 
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backwards  shifts  of  center  of  gravity.  The  problems  of  poor 
response  characteristics  and  controls  going  out  of  bounds 
only  occured  at  flight  conditions  well  removed  from  the 
nominal  and  posed  no  problem  while  near  the  nominal  flight 
conditions.  Furthermore  these  problems  could  be  alleviated 
by  using  a different  feedback  matrix  or  by  incorporating  the 
control  bounds  in  the  controller  itself. 

Probably  the  most  serious  deficiencies  to  this  method 
is  the  unfavorable  response  characteristics  at  flight 
conditions  away  from  the  nominal.  The  technique  most  widely 
used,  currently,  to  bypass  this  problem  is  to  have  many 
precalculated  feedback  gain  matrices  available,  and  to  just 
use  one  in  accordance  with  the  flight  condition  the  matrix 
most  closely  corresponds  to.  This  is  a practical  solution 
but  seems  somewhat  awkward.  A more  recent,  and  perhaps  more 
promising  concept  involves  the  use  of  adaptive  control. 
This  method  automatically  adjusts  the  gain  matrices  in 
accordance  with  the  changing  flight  conditions.  This  type 
of  control  has  been  applied  to  vertical  takeoff  and  landing 
(VTOL)  aircraft  with  favorable  results  and  quite  possibly  is 
the  method  of  the  future. 
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APPENDIX  A 

DESIGN  OF  A LOW  ORDER  OBSERVER 


Given  a system  of  the  form: 

x = Ax  + Bu  A[n  x n] 


B [ n x m ] 


y = Cx 


CLp  x n] 


where  (A,C)  is  observable  and  C is  of  full  rank  ( Ref.  4). 
Choose  a transformation  matrix  T = [C  M]  where  M is  an 
(n-p)  x n matrix,  such  that  T is  nonsingular  (often  chosen 
so  that  T simply  interchanges  the  order  of  the  state 


variables 


transf ormation : 


( i . e . 


:2»x1»x4,...]) 


Using  this 


z = fz i I = Tx  = [ y 
z2  Z 

z = TAT"^z  + TBu 


A 21  A22 


= TB  = B. 


z-|  - A 1 1 z i + ^ ] 2^2  + B^u 
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2)  Z2  = A2 i z i + ^22z2  + B2U 

Since  z-j  = y is  already  known,  we  only  need  an 
observer  for  Z2.  Now  define: 

3)  £ 2 = ^22z2  + ^21z1  + + K C A 1 2 2 2 ~ A 1 2Z2  ) 

and  call  e = §2  - Z2 

4)  = $T(z2  ~ z2)  ' (A22  " KA12)(z2  - Z2) 

The  eigenvalues  of  this  can  be  placed  such  that  the 
error,  e = 0 in  the  steady  state.  This  is  possible  since 
(A,C)  is  observable  which  implies  (A22.A12)  is  controllable. 

From  1)  and  3)  we  see  that: 

2 2 = A 2 2 z 2 + ( A2 1 z 1 + B 2 u ) + K ( z 1 - A 1 1 z 1 
— B 1 u — A 1 2Z2 ) + ^22  — KA 1 2 ) (Kz  1 ) 

•^($2  ~ Kz-j)  = (A22  ” KA 1 2 H ^2  “ Kz  1 ) + (A21  - tc  3 1 1 
+ (A22  - KA12)K  )zj  + (B2  - KB-j  )u 


then 


where : 


w 2 = z 2 ~ Kz 1 


W2  = G 1 w2  + G2y  + G^u 


G 1 = A22  " KA12 

G2  = A 2 1 “ KA 1 1 + G 1 K 

G3  = B2  - KB  1 


This  yields  an  observer  system  as  shown  in  Fig. 9,  and  x 
may  now  be  obtained  from  the  equation: 


x = T" 


'N 
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APPENDIX  B 
MINIMUM  PRINCIPLE 


Given  a system  of  the  form  : 
x = Ax  + bu 
y = Cx 

and  a performance  index: 

/OO 

( xTQx  + uTRu ) dt 

u 

Where  Q and  R are  positive  semidefinite  and  positive 
definite,  symmetric  cost  matrices.  Using  the  techniques  as 
discussed  in  Athans  and  Falb  (Ref.  1),  form  the  Hamiltonian: 
H = L + pTx  = (1/2)(  xTQx  + uTRu  ) + dT(Ax  + Bu) 

From  the  Minimum  Principle,  for  unconstrained  u 
(hopefully  u will  stay  within  constraints): 

0 = fu  = uTr  + PTB 
u = R“1BTp 

P = “|x  = -Qx  ~ APT 
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If  we  now  assume  p = Kx  (as  discussed  in  Athans  and 
Falb)  where  K is  time  invariant  since  A,B,Q,  and  R are  all 
time  invariant,  we  find: 
u = R_13tKx 

From  the  costate, p,  equation , then 

KA  + ATK  + KBR_1BTK  + Q r 0 
Which  can  be  solved  for  a unique  K,so 
u = Fx  where  F=R“13Tp 
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APPENDIX  C 

COMPUTER  PROGRAMMING 


To  get  analytical  results  for  all  of  the  equations  just 
derived  computer  programs  had  to  be  written  to  perform  all 
of  the  integrations  and  other  related  operations  needed. 
For  reasons  of  time  and  accuracy,  all  preliminary  simulation 
was  done  on  the  DEC  10  digital  computer.  Final  results  were 
then  tested  on  the  hybrid  computer,  where  the  AD-5  analog 
computer  represented  the  real  time  plant  (an  airplane)  which 
was  being  controlled  by  the  digital  PDP-11  digital  computer. 


ft  | 

\ ■ 


1 .DEC-10 

All  control  testing  and  sensitivity  analysis  was  done 
on  the  DEC-10.  Because  of  the  size  of  the  program  and  the 
need  for  versatility  of  input  data,  an  interactive  format 
was  utilized.  This  method,  of  having  the  operator  respond 
to  different  options  (e.g.  initial  conditions ) .helped 


facilitate  debugging  of  the  program  also.  Furthermore  this 
method  made  it  possible  to  study  any  flight  condition  or 


— 


possibly  any  aircraft  conf iguration  within  reason,  by  a 
simple  response  to  a parameter  change  option.  The  only  true 
shortcoming  involved  here  was  that  the  program  did  not  have 
the  option  of  generating  feedback  matrices  (this  was  done  on 
Linsys,  using  the  Riccati  package)  so  the  responses  to 
different  conditions  (other  than  the  initially  chosen  one) 
were  suboptimal  in  some  sense. 

All  of  the  interactive  programming  and  condition 
organization  was  done  with  one  very  large  main  program, 
ACUTI.  This  program  would  ask  for  the  desired  flight 
conditions  or  the  desired  data  analysis  (such  as  sensitivity 


or  observer  analysis) 

and 

would 

then  make 

calls  to 

the 

various  subprograms 

needed 

to 

facilitate 

these. 

The 

subprograms  would  then 

execute  the 

dif  f erent 

commands 

such 

as  for  integrations  or  plots. 


2. HYBRID 

The  AD-5  analog  computer  had  the  nonlinear  aircraft 
plant  equations  patched  into  it,  thus  simulating  a real  time 
airplane.  This  required  a lot  of  manipulation  and  scaling 
due  to  the  limited  amount  of  hardware  available,  and  due  to 
saturation  restrictions.  To  help  set  up  and  test  this, 
several  PDP-11  programs  were  used.  Again,  here,  the 
programs  were  set  up  interactively,  so  theoretically,  any 


Li 


flight  conditions  could  be  simulated.  But  here,  again,  due 
to  scaling  and  hardware  limitations,  there  was  actually  only 
a limited  range  of  variations  possible.  For  accuracy  and 
speed  of  setting  up,  a subroutine,  POTREV,  was  used  to 
calculate  and  set  all  pot  values,  automatically,  according 
to  what  parameters  were  desired.  The  analog  diagram  is 
shown  in  Fig . 1 0 . 
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